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Abstract 

One-particle inclusive CP asymmetries in the decays of the type B 

are considered in the framework of a QCD based method to 
calculate the rates for one-particle inclusive decays. 



PACS numbers: 11.30.Er, 13.25.Hw 



1 Introduction 



One of the main goals in B physics is a detailed study of flavor mixing, 
which is encoded in the Cabibbo-Kobayashi-Maskawa (CKM) matrix of the 
standard model. In particular, the violation of the CP symmetry, which 
the standard model describes by a nontrivial phase in the CKM matrix or 
equivalently by the angles of the unitarity triangle, will be investigated. 

Typically CP asymmetries are expected to be large in some of the ex- 
clusive nonleptonic B decays which, however, have only small branching ra- 
tios. Examples are the determination of f3 from B J/ip Ks and of a from 
B ^ 7171. In addition, in these exclusive nonleptonic decays it is very hard 
to obtain a good theoretical control over the hadronic uncertainties, in par- 
ticular due to the presence of strong phases. 

On the other hand, inclusive decays have large branching fractions but 
typically smaller CP asymmetries than exclusive decays One may use 
parton hadron duality to obtain a good theoretical description. This has been 
studied by Beneke, Buchalla and Dunietz who set up a theoretically clean 
method to calculate the CP asymmetries in inclusive B decays 0. They 
still find sizable CP asymmetries, but their measurement would require to 
identify charmless final states inclusively, which is not an easy task. 

One-particle inclusive decays lie somehow between these two cases. This 
class of decays still has large branching fractions and some of the expected 
CP asymmetries are sizable. Furthermore, a measurement of these decays is 
feasible. 

For one-particle inclusive decays of the type B —>■ D^*^X, a QCD based de- 
scription has been developed recently, exploiting factorization and the heavy 
mass limit for both the b and the c quark ||^. Since the expansion parame- 
ters are Aqcd/("^6 — ^c), and ^^(mc), corrections to the leading term 
could be fairly large, in the worst case of the order of 30%. Using this 
method, which unfortunately is not completely model independent, we com- 
pute mixing induced time-dependent and time-integrated CP asymmetries 
in the framework of the standard model. 

In view of the considerable uncertainties due to an unknown strong phase, 
our method cannot yet be used for a competitive determination of the CP 
violation parameters, in particular compared to a measurement of sin (2/5) 
in the "gold-plated" channel B —>■ J/ipKg. However, it can be used as an 
estimate of the one-particle inclusive CP asymmetries, for which we shall 
use present central values of the CP angles /? and 7 |Q. Compared to fully 
inclusive methods, the advantage is that we can predict asymmetries for the 
various spins and charges of the ground-state charmed mesons separately. 
This is certainly a worthwhile task, in particular since we are not aware of 
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any previous prediction for these asymmetries, not even in the context of 
quark models. 

After introducing our notations for B mixing in Sec. we calculate the 
relevant matrix elements in Sec. ^ and model the form factors in Sec. The 
numerical results are given in Sec. |. 



2 CP asymmetries in B — > 

In Wigner Weisskopf approximation the time evolution of an initially pure 
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where AM = Mh — > and Ar = r^^ — r^ < are the mass and 
width differences between the mass eigenstates \Bh> = p\B^> + q\B > and 

\Bl> = p\B°> -q\B°>. 

The quantity q/p is given in terms of the off-diagonal elements of the 
Hamiltonian H = M — iT /2 of the neutral B meson system 
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In fact, T12/M12 = 0{ml/m^) is very small and hence q/p is to a good 
approximation a phase factor. 

The time-dependent rate for the decay of a i? meson into a set of final 
states I/) = J2i \fi) can be written as 

r[B{t)^f] = Jd<l>d27r)'6'{pB-Pf^){B{t) \H,s\ f^) {f.\H,s\ B{t)) 
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where d(j)i is the phase space element of the state \fi) and 

Ilf = T.J (/^l 
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is the projector on the set of final states. Note that both an exclusive final 
state as well as inclusive states can be treated in this way. Even differ- 
ential distributions can be considered if the phase spaces d(f)i are not fully 
integrated. 

The CP asymmetries we are going to consider are of the type 
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which involves the CP conjugate set |/> of final states. 

Up to here the discussion is completely general. In the following we shall 
use the above formalism to compute the CP asymmetries for one-particle 
inclusive final states, for which the projector reads 



= ^ \XY) {XY\ , 
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where Y can be a D or a D meson. Since the sum runs over all possible 
states X, the CP conjugate of the projector is 



X 



Inserting the time-dependent states (|I|) we obtain 
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where the matrix elements are defined by 
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The AB = 2 transition matrix elements representing the interference between 
the mixed and the unmixed amplitudes are related by CPT symmetry, such 
that 
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The direct CP asymmetries in these processes are expected to be tiny. 
In fact, using the method described in Ref. 0, they turn out to be of higher 
order in the 1/m expansion. Hence we have 

Ty := = = T{B ^ YX), Ty := if^, (12) 

Ty = Ty. (13) 

Inserting the time-dependent decay rates in Eq. (H) and neglecting both 
the width difference and a, such that q/p becomes a phase factor, we obtain 
for the time-dependent CP asymmetries 

sin (AMt) Im {jTy} 

from which we get the time-integrated asymmetry 

2a;Im{2Ty} 

(2 + X^) Ty + X'^ Ty 

where x = AM/T is measured to be x = 0.73 ||^. 



3 Transition matrix elements 

In order to compute the CP asymmetries, one has to evaluate the matrix 
elements in Eq. (0). The total rates Ty have already been discussed in 
Ref. [0, so we only need to calculate the interference term Ty. 

The relevant pieces of the effective Hamiltonian contributing to this inter- 
ference are (ub)v-A{dc)v-A and (cb)v-A{du)v-A interfering with each other 
and {cb)v-A{dc)v-A interfering with itself, so Ty is a sum of the two contri- 
butions 
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Fierzing the operators into the form {db)v-A^c)v-Ai {,db)v-A{cu)v-A 
and {db)v-A^c)v-A one can reproduce the inclusive results of Ref. 0. In 
order to evaluate the interference term for the one-particle inclusive case, we 
use the method developed in Ref. H. It is based on factorization, which 
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holds to leading order in the I/Nq expansion, where Nq is the number of 
QCD colors. Thus we can write the interference terms as products of two 
tensors 
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mi'{^-l^)c)\D^\pD)X') (D«(p^)X'|(c7^(l-75)g)|0). 



B\pb) (^7^(1-75)6) 



^^(PD, Q) = E(2^)'^'(Q -PD- PX') 



The tensor K^y{pBi Q) is fully inclusive and one can perform a standard 
short distance expansion. The resulting AB = 2 matrix element can be 
parameterized by the decay constant Jb of the B meson and the bag factors 
B and Bs for the axial vector and the scalar current, respectively. 

The other tensor P^^IpbjQ) involves a projection on a one-particle in- 
clusive charmed meson state and hence we cannot perform a short distance 
expansion. We proceed along the same lines as in Ref. [Q, where the rates 
for wrong charm decays have been modeled. Heavy quark symmetry yields 
the Dirac matrix structure 



Pr{PD,Q) oc H^,.,{pn)r{l-l,)®Y{l-l6)H^i*){pD) 
where the representation matrices for the charmed mesons are 
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In principle, all possible contractions of the light quark indices may con- 
tribute, giving rise to several form factors. For a first estimate, it is sufficient 
to use only the simplest one of these contractions. 



P'^{pr,^Q)=2i^b([Q-pu) 



1 2 



Tr{|^z5 7^(l-75) (^-I^d) 7^(1-75)} 

(23) 

corresponding to a replacement of the D'^*^X final state by a pair of free 
quarks, rescaled by an operator- and decay-channel-specific form factor /^y, 
where Y is one of the ground state D mesons. In the following, we call this 
contraction "partonic." 
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Using this ansatz and the heavy mass hmit, the transition matrix elements 
read 



= -^^^^{V,,V:,Y \C,\' VT^z[{l-Az)B + 2{l + 2z)Bs]ly, (24) 

= -^^^^v;,K*.K.K*.|c^i|'(i-^)'^ (25) 

where z = {mc/mi,Y and Ci is the Wilson coefficient of the effective Hamil- 
tonian in the notation of Ref . . Equations (|^) and (|25|) correspond to the 
expression for the width difference of neutral heavy meson systems [^]. 

In the standard CKM parametrization, the phases of the transition matrix 
elements are 

(26) 
(27) 
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such that 



Im I ^^5- 1 = sin(2/3) |r,| + sin(2/? + 7) |T„| . (29) 



4 Modeling the form factors 

We assume that the form factors fqy do not vary strongly over the accessi- 
ble phase space and hence we approximate them by constants. For the case 
q = c, these constants have been fitted to the wrong charm yield in B de- 
cays . Operators analogous to the case q = u are Cabibbo suppressed when 
calculating wrong charm rates, so they did not appear in Ref. . Assuming 
that all charm quarks eventually hadronize to D mesons, we use 

fuDO + fuD+ = 1. (30) 

To resolve the spin and charge counting, we first discuss the heavy mass 
limit where the pseudoscalar and vector charmed mesons form a degenerate 
ground state doublet. The decay of vector to pseudoscalar mesons will be 
discussed below. In the following, Ddh refers to those D mesons that do not 
result from D* decays, and D^*^ can be either Ddir or D*. 

As long as the light quark spin indices of the D^*'^ meson representation 
matrices are contracted with each other, Eq. ( PT| ) reproduces the naive spin 
counting 

fqD'O = '^fqDO.^^ fqD*+ = '^fqD+^- (^1) 
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Figure 2: Topology yielding joal^ Figure 3: Interference topology 
in Eq. (H). for Eq. 

Different contractions yield results of comparable size. The experimental spin 
counting factor appears to be smaller by roughly a factor of two [Q]. Since 
this effect is not yet understood, we treat it as an uncertainty. 

Concerning charge counting, we argued by isospin symmetry that in 
the case g = c we have 

fcDMO = fcDM+- (32) 

In the case q = u, two topologies can contribute to the decay amplitude: the 
charm quark can either hadronize with the u quark from the weak effective 
current, in which case the isospin of the state |X) is = 0, or with a u 
or d quark from vacuum, which contains both Ix = and Jx = 1 contribu- 
tions. In the case Ix = 0, both amplitudes can interfere, so there are three 
contributions to the decay rate 

fuD(*)o = \ai + = \aif + \a2f + 2Re{a*a2} ^gg^ 

fuD(*)+ = , 

see Figs. 

One might doubt whether using the partonic contraction given in Eq. ( P5| ) 
is justified for all the topologies, as it appears to correspond to the topology 
in Fig. H, while the topology in Fig. |I| should rather be described by the 
contraction 

P!;''{PD.Q) oc Tr{:ff^„(p^)7^(l-75)}Tr{7''(l-75)i^Dw(PD)}. (34) 
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This is not a problem for three reasons. First, we do not claim to be 
able to accurately model the matrix element, but we only give the simplest 
possible ansatz by rescaling the partonic result. In particular, it is clearly 
not yet feasible to model particular contributions individually. We only use 
the three topologies to estimate the integrated relative magnitudes of the 
two main contributions and to bound the magnitude of their interference 
term. Secondly, neither the time- dependent nor the time-integrated asym- 
metries depend on the choice of the contraction unless studied differentially 
in the momentum of the charmed meson, which so far we do not attempt to 
do. Finally, as noted in Ref. the choice of the wrong charm contraction 
appeared to have little influence even on differential observables. 

The topologies in Figs. |I] and ^ also occur in wrong charm production 
in B decays. Figure D corresponds to the process B Dl*^^X, Fig. || to 
the process B D^*'>X, where D^*'> can be either or Both 

contributions are experimentally known to be of similar size, i.e., (10±2.^ 
§ and (7.9 ± 2.2)% 0], respectively, such that 



ai 



2|a2| 



(35) 



The relative phase of the two contributions is unknown. Therefore, although 
it may be large, we have to treat the interference part as a theoretical un- 
certainty. This is acceptable since the q = u contribution is smaller than the 
q = c contribution according to 
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Off the heavy mass limit, D* D decay has to be taken into account. 
In the same way as in Ref. [H], we get 



f,D+ = hot, + Br {d*+ ^ D'-X) f,n. 
f,DO = hnl^ + /,D*o + Br {d*+ ^ D^X) /, 



(37) 



qD* 



The coefficients obtained from Eqs. ([30|)-(|37D and Ref. are summarized 
in Table |^. The ranges given result from varying the spin counting factor 
in Eq. ([3l| ) from 3 down to 3/2 and the interference in Eq. (|33|) from the 
central value of vanishing interference to full constructive and destructive 
interference. 
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^^\operator 


q = 


c 


Q 


= u 


channel 


central 


r. to 


central 


range 


^ dir 


2/16 


0.2 


1/16 


0.04-0.34 




2/16 


0.2 


3/16 


0.34-0.04 


fqD''+ 


6/16 


0.3 


3/16 


0.09-0.64 


fqD'O 


6/16 


0.3 


9/16 


0.64-0.09 


fqD+ 


4/16 


0.3 


2/16 


0.07-0.51 


fqDO 


12/16 


0.7 


14/16 


0.93-0.49 



Table 1: Operator- and channel-specific form factors. 



5 Results 

We have computed the parameters for the time-dependent CP asymmetries 
as well as the time-integrated asymmetries. We have inserted recent val- 
ues for sin 2/? = 0.75 and 7 = 68° [Q. In addition, we use Vcb = 0.04, 
Vub = 0.08 Vcb, z = 0.09, X = 0.73, fe = 180 MeV, Br {D*+ D^Y) = 
1 - Br D+Y) = 0.683 and Ci = B = Bs = I. The results of the 

calculations can be found in Fig. ^ and Table ^. 

To assess the uncertainties involved in Fig. |^, note that according to 
Eq. ( [T^ ) the shapes of the time-dependent asymmetries are determined by 
the ratios of the wrong to right charm rates Fy/Fy. We checked numeri- 
cally that the shapes would hardly change even if these ratios were off by 
30%. The dominant contribution to the uncertainty of the amplitudes arises 
from the transition matrix elements Ty and is directly proportional to the 
uncertainties of the time-integrated asymmetries given in Table 0. 

Suppose perfectly tagged 5° decays are recorded in an experiment. In 
order to establish the asymmetry in a channel with a branching ratio b on 
the 3a level, 

- > = -L= (38) 

has to be satisfied. The necessary numbers of tagged B^ decays are given 
in the last column of Table |^. Since the asymmetry tends to be roughly 
inversely proportional to the branching ratio by Eq. (|15|), we obtain from 
Eq. (HD 

AT cx (X 6, (39) 

such that rare channels are advantageous for observing one-particle inclusive 
asymmetries. 
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Figure 4: Time- dependent CP asymmetries in 5° — DX for pseudoscalar 
(above), vector (below), charged (left), neutral (right), right charm (solid), 
and wrong charm (dashed) jD mesons. 
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decay 


Br i 


A 


A range 


necessary 


channel 


(%) 


(%) 


(%) 


i?o decays 


5° D'X 


29.1 


0.16 


0.15-0.29 


6.000.000 




31.8 


0.58 


0.59-0.46 


400.000 


fiO ^ D+X 


2.2 


0.58 


0.54-1.04 


6.000.000 




5.7 


1.53 


1.56-1.23 


350.000 




46.8 


0.16 


0.12-0.23 


4.000.000 




(0) 


(20) 


(21)-(10) 


> 80.000 


5° ^ D*+X 


2.5 


0.61 


0.45-0.89 


5.000.000 




2.5 


4.17 


4.40-2.20 


100.000 



Table 2: Branching ratios, integrated CP asymmetries and numbers of nec- 
essary tagged decays for the one-particle inclusive B^ D^*^X decay 
channels. Concerning D , see the text. 

The channel B^ D*^X deserves a further comment. Looking at Fig. ^, 
there is an obvious problem at small proper decay times. The reason for this 
problem is that we have discussed all the rates only to leading order in the 
combined 1/Nc and I/tuq expansions. However, this leading term vanishes 
for the channel D*^X and thus subleading terms become relevant. 

On the other hand, the numerator Ty of the CP asymmetries is given by a 
matrix element of a dimension six operator and hence is suppressed compared 
to the leading terms of most of the rates. In other words, while in most of 
the rates the asymmetries are of subleading order Jb/^'bi this is not the case 
for the channel B^ D*^X. 

Unfortunately we cannot compute this possibly large asymmetry, since 
this would involve to compute subleading terms for the decay rate. Hence we 
try to estimate the asymmetry by varying Br (5° ^ D*^X) in Eq. (||) and 
show the reaction of the asymmetry in Fig. |^ and of the necessary number 
of tagged B^ events in Fig. ^. The wrong charm asymmetry is practically 
unaffected by Br (^B^ D* X^ since the pole occurs near four average life- 
times where most of the B mesons have already decayed, but the right charm 
asymmetry turns out to be extremely sensitive. Therefore we cannot predict 
the latter quantitatively, but it can be as large as several percent, and it will 
be measurable with a few 100 000 tagged B^ events. 
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Figure 5: Time- integrated asym- 
metry in D*^X as a func- 
tion of Br (5° ^ D*°X). 



Figure 6: Necessary number of 
tagged E° events in B° D*^X 
as a function of Br(i?° ^ Z}**^X). 



6 Conclusion 

Motivated by the work on fully inclusive CP asymmetries and the question 
how to measure them, we studied one-particle inclusive CP asymmetries. In 
the final state only a D*-*-* meson has to be identified and thus they are ex- 
perimentally more easily accessible than the fully inclusive CP asymmetries. 

We have used a similar method as in in Ref. to calculate the time- 
dependent and time-integrated CP asymmetries for one-particle inclusive 
B D^*^X decays. It turns out that, as in Ref. 0, one cannot avoid to 
introduce some model dependence. Furthermore, there is also some depen- 
dence on an unknown relative phase, which we treat as an uncertainty. Due to 
these uncertainties we cannot expect our method to compete with proposed 
methods using "gold-plated" channels for determining CKM parameters, but 
we can still give estimates for the expected CP asymmetries of the different 
ground state D mesons. 

For most of the asymmetries we find results of a few 10"'^, but some are 
expected to be as large as several percent. These effects should be observable 
at the B factories. The channels involving right and wrong charm neutral 
vector mesons turn out to be most promising: they are expected to have the 
largest asymmetries, and the theoretical method yields the best results for 
the production rates and spectra of the vector mesons 0. 
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